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We consider energy absorption in an externally driven complex system of noninteracting fermions with 
the chaotic underlying dynamics described by the unitary random matrices. In the absence of quantum inter- 
ference the energy absorption rate W(t) can be calculated with the help of the linear-response Kubo formula. 
We calculate the leading two-loop interference correction to the semiclassical absorption rate for an arbitrary 
time dependence of the external perturbation. Based on the results for periodic perturbations, we make a con- 
jecture that the dynamics of the periodically-driven random matrices can be mapped onto the one-dimensional 
Anderson model. We predict that in the regime of strong dynamic localization W(t) oc h\{t)/t 2 rather than 
decays exponentially. 
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1. Introduction. Last years had revealed an in- 
creasing interest [fl [21 El ED to the time-dependent ran- 
dom matrices, arising from the field of condensed matter 
physics. The natural way to study a complex quantum 
system is to couple it to an external field ip which enters 
the Hamiltonian H[(p\ = Ho + Vip as a parameter and 
can be controlled at will. Applying a time-dependent 
perturbation ip(t) gives access to quantum dynamics 
of the many-electron wave function governed by the 
Schrodinger equation i d^(t)/dt = H[ip(t)] If the 

perturbation frequency and the relevant energies (e.g., 
the electron temperature) are smaller than the Thou- 
less energy in the sample then it is possible to apply 
a universal description in terms of the random-matrix 
theory (RMT) of an appropriate symmetry 0. The re- 
sulting time-dependent theory is specified by two model- 
dependent quantities, which should be determined mi- 
croscopically H3 : the mean level spacing A and the sen- 
sitivity of the parametric spectrum ((dEi/ 'dip) 2 ) to the 
variation of the control parameter ip. 

The crucial quantity characterizing quantum dy- 
namics of the system is the energy absorption rate 



W(t) 



d(E(t)) 
dt 



(1) 



and its dependence on the form of the external pertur- 
bation (p(t). [In Eq. (E(t)) is the expectation value 
of the total energy of the system.] The standard ap- 
proach to calculation of W is based on the Kubo linear 
response theory which expresses the energy absorption 
rate in terms of the matrix elements of dH/dt. For the 



standard Wigner-Dyson random matrix ensembles one 
finds 

Wo = ^C f Mv\ (2) 
where v = dip/dt is the perturbation velocity, 



Cp(0) 



1 

A 2 



dEi\ 



/3A 2 



(3) 
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is the level velocity autocorrelation function, with Ei[ip] 
being the adiabatic levels of an instantaneous Hamilto- 
nian, and (3 = 1 or 2 for the orthogonal (GOE) or uni- 
tary (GUE) symmetry classes, respectively. The Kubo 
dissipation rate @ is ohmic as it scales ex v 2 regardless 
of the system's symmetry. 

The semiclassical result (j5J| was obtained neglect- 
ing quantum phenomena in dynamics. There are two 
types of interference effects which may invalidate the 
semiclassical description. The first one is related to the 
condition of continuous spectrum implicitly assumed in 
evaluating the Kubo commutator. For a closed system 
the Kubo formula J2J) can be applied only at sufficiently 
large v ^> vk ~ A 2 / v/(7g(0) when the spectrum is 
smeared by nonstationary effects. For small v <C vk 
the dynamics is adiabatic and dissipation is due to rare 
Landau-Zener transitions between the neighboring lev- 
els. In this case the energy absorption rate becomes 
statistics-dependent with W ~ u^/ 2+1 . The second 
interference effect comes into play for re-entrant per- 
turbations when the system is being swept through the 
same realization of disorder many times. For a certain 
type of time-dependent perturbations, destructive inter- 
ference in the energy space may lead to dynamic local- 
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ization jEJ] and hence to the vanishing of the absorption 
rate. 

Recently the first quantum interference correction to 
the Kubo dissipation rate (J2J) for the orthogonal symme- 
try class was considered, taking into account both the 
original discreteness of the spectrum [3J and the effect 
of weak dynamic localization 0]. The one-loop rela- 
tive correction to Wq contains a dynamic cooperon and 
evaluates either to a positive number ~ {v/vk) 2 ^ 3 for a 
linear bias ip — vt [3] or to a negative and growing in 
time correction oc —y/i for a monochromatic perturba- 
tion switched on at t = @] (in this case the dynamic 
localization effect is the most pronounced). 

The purpose of the this Letter is to study the quan- 
tum interference correction to Wq for the unitary sym- 
metry class, that involves evaluation of the two-loop di- 
agrams made of dynamic diffusons. We will derive the 
general expression for 8W(t) [Eq. (|21[l ] valid for an ar- 
bitrary time dependence of p(t) and then discuss the 
limits of linear and (multi-) periodic perturbations. 

2. Description of the formalism. Quantum dy- 
namics of time-dependent unitary random matrices can 
be conveniently described by the nonlinear Keldysh a- 
model derived in Ref . ■ The effective action (with the 
weight e~ 5 ) 



S[Q] = -TtEQ- 



7T 2 C M (0) 



Tr[^,Q] 



(4) 



is a functional of the Q field acting in the Keldysh (Pauli 
matrices <7j) and time spaces. In Eq. (@J the operators 
E and p have the matrix elements Ew = iStt'df and 
<ptt' = Stt'p{t r ), and C u (0) is the level velocity autocor- 
relation function defined by Eq. with = 2. 
The saddle point of the action (@J is given by 



htt> = 



. 



(5) 



with the distribution function F^ satisfying the kinetic 
equation 

(d t + ft,) f£) = -r Mt) - p{t')Y f$\ (6) 

where we denoted T — irC u (0)A. 

The whole manifold of the Q matrices can be 
parametrized as 



= Up 1 PU F , P=U- 1 a 3 U, 



(7) 



where the matrices U are unitary, so that P is a Hermi- 
tian field, whereas all non-Hermiticity is located in the 
matrices 



(U F )t 



, -Ste 



(8) 




(b) 




(c) 




Fig. Two-loop diagrams for the distribution function 
F, corresponding to the terms of Eq. 1151 . Solid lines 
denote the diffusons. 

[in particular, the standard saddle point (J5J corresponds 

t0P = CT 3 ]. 

For perturbative calculations we choose the standard 
rational parameterization of the P matrix, 



P = <r 3 (l + V/2)(l-V/2)-\ 



(9) 



which has the unit Jacobian dP/dV = 1. The matrix 
V anticommuting with a 3 is given explicitly by 



V = 



d\ 
-<ft ' 



(10) 



with the matrix d acting in the time space only. Its bare 
correlator inferred from the Gaussian part of the action 
has the form: 

2A 

(d t+t _(% +t ,_)o = —5(r 1 -7 1 , )V tl (t,t'), (11) 

where we have denoted t± = t ± ry/2, t'± = t'±rj' /2, and 
introduced the free diffuson propagator 1101 

V r] {t,t') = 9{t-t') expj-^ T[<p(T+) - p{t-)] 2 dr\ . 

(12) 

Physical quantities are contained in the average 
(Q) = J Q e~ s ^ DQ. Due to causality, (Q t t>) shares 
the structure of the Eq. (JSJ but with the saddle-point 
distribution F^ substituted by the exact distribution 
F. The energy absorption rate can be calculated as 0] 



W(t) 



limd t d n F t+n /2^- n /2- (13) 



A >7->o 



3. Perturbation theory. Expanding the Keldysh 
(upper-right) block of the matrix Q in terms of the dif- 
fusons d with the help of Eqs. (f^Zf) — (| 1 0|) one obtains the 
perturbative series: 



F 



p(o) _ <d> _ (F^dU + dd^F^) + (dd^d) 



(14) 

The two-loop correction to the distribution function is 
given by three pairings: 

°F= n o + 3 > ( 15 ) 
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shown diagrammatically in Fig. ^ The other possi- 
ble pairings vanish due to causality of the theory. In 
Eq. JTSJl, the vertices S( m > come from expansion of the 
action (@J to the order W rn . In the rational parametriza- 
tion they are given by the following expressions: 



7lT 

2A 



J V 1 2 </?34-FisP ^32 ^34^14 > 



(16) 



S {5) = ^ ^12(^4 + ^56)^12^23^43^45^65^61 



(17) 



8A 



(<9 5 + d 6 )d 56 dj 6 d 78 d* 58 



16A 



(v>56 + ^58 + ^67 + ^78 



^57 _ ^68 ) d 5 6d%d 78 d : 



58 



(18) 



The the terms not included in Eqs. I|17l) and (|18fl do not 
contribute to the pairings shown in Fig. ^ In writing 
Eqs. JTSJ-lQSJ we use d the concise notations Fy = F^u , 
dij = dtitj j and tpij = ip(ti) — <p(tj), with integration be- 
ing performed over all time arguments involved. 

The diagrams (a) and (b) shown in Fig. ^ contain a 
loose diffuson 2 which couple d to the rest of the dia- 
gram. As a result, the corresponding correction to the 
distribution function i 7 t+^/2,t- J) /2 can be written as 



SF 



(ob) 



t+ri/2,t—ri/2 



J dt'dt" V„(t, t') E(t', t", v ) F t „ +v/2 ^_ v/2 , 



(19) 



where t' is the "center of mass" time at the right end 
of the loose diffuson, and E(t',t",n) is a complicated 
expression denoting the rest of the diagram. The corre- 
sponding correction to the energy absorption rate given 
by Eq. (|13|l simplifies to 



5W^(t) = -i Urn iJLl f dt"S(tX,v), 
A 17— >Q or/ n J 



(20) 



where we employed Eq. (|12[) and used the asymptotics 
F t+t _ ~ l/(wr??) at rj -> 0. 

Contrary, the diagram (c) in Fig. ^ does not contain 
a loose diffuson and cannot be represented in the form 
(|2U[I with already taken derivative with respect to the 
external time t. 

Finally, it is worth mentioning that the diagram (a) 
is completely canceled against the part of the diagram 
(b) which contains the time derivative originating from 
the first term in Eq. ((TBI) . 



As a result of straightforward but rather lengthy cal- 
culation one ends up with the general expression for 
the two-loop correction to the Kubo dissipation rate 10) 
valid for an arbitrary tp(t): 

TA S 1 f°° 
8W(t) = -—7 lim - — / dxdydz 
y ' 2tt 2 »j-o dijij J y 

X ( - 2r^5 6 937 8 J V?i2</?34 



X V. 



r)+x+y 



XV XV 
t -, t --Z 

2 2 2 2 



x V. 



X z 



X z 



\ { 2 2' ' 2 2 V 



, y z y z 

X TJr)+y — z I t — — — — , t 



x), (21) 



where t\^ = t± — x — y — z, t% = t + — z, t± = t- — y, 
£5 = t + — x — z, te t 7 = tip, and tg = t- — x — y. In 
Eq. jSJJ) the term with d/dt describes the contribution 
of the diagram (c) while the rest is the contribution of 
the diagrams (a) and (b) . Thought the derivatives with 
respect to rj and t can be easily calculated with the help 
of Eq. I|12l) we leave them unevaluated in order to keep 
the simplest form of the expression. 

4. Linear case. We start the analysis of the 
general formula (|21|l with the case of a linear bias 
ip(t) = vt. Then the dynamic diffuson (|12|l is given 
by T> v {t u t 2 ) = 0(*i - * 2 )exp{-nV(ii - fc)} where 
n = iTv 2 ) 1 / 3, is the dephasing rate due to the time- 
dependent perturbation 3 . Since the diffuson T> n (t\ , t 2 ) 
depends only on t\ — fa, the integrand in Eq. l|2"T|> does 
not depend on t and the corresponding time derivative 
describing the contribution of the diagram (c) vanishes. 

The product of three diffusons in Eq. I|21(l is an even 
function of rj, hence ^-dependence should be taken into 
account only in the terms ipij . The resulting expression 
becomes 



SW = 



n 6 A 



dx dy dz (—x 2 + 5xy) 



xexp{-r! 3 (x + y)(y + z)(z + a; )}, (22) 



where we employed the symmetry between the integra- 
tion variables to simplify the final expression. 
The integrals in Eq. I|22|) are given by 



dx dy dz 



■ T V. 



-(x+y)(y+z)(z+x) _ 



r 2 (i/3) 

48 : 



(23) 

leading to a surprising cancelation of the two-loop quan- 
tum correction in the unitary case mentioned in Ref. [3] . 

It is also instructive to consider the case of the linear 
perturbation switched on at t — 0: ip(t) = 9(t) vt. Here 
the term with d/dt in Eq. i|21|) is generally nonzero but 
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it is small in the most interesting limit fit 3> 1. The 
time-dependent SW(t) is then given by Eq. (|22|l where 
the region of integration is now bounded from above by 
the condition x + y + z < t. The correction to the total 
absorbed energy becomes 

tt 6 A f°° 

SE(t) = — 5— / dx dy dz minfizi + y + z, t) 
* Jo 

x {-x 2 + 5xy) e -n 3 (*+y)(y+z)(z+x) ^ 

The integrals with x 2 y and xyz converge while the in- 
tegral with a: 3 diverges logarithmically. Therefore, at 
fit > 1 

5E(t)~-^bi(nt). (25) 

Thus, the two-loop quantum correction, though vanish- 
ing for a linear perturbation, leads to a long-time mem- 
ory effects near the points of discontinuity of dtp/dt. 

5. Periodic case. Now we turn to the case of peri- 
odic perturbations switched on at t = 0. To simplify cal- 
culations we will consider first the simplest example of a 
monochromatic perturbation, (pit) = sinwi. Then 
the dynamic diffuson 1)12(1 acquires the form: 



D v (t,t') =0(i-i')exp 



-2rsin 2 — 



, smujft — t 1 ) , ,. 
t — t A cosw(i + t 1 ) 



(26) 



It is convenient to calculate the two contributions to 
Eq. (EB, SW^(t) and 5W^(t), separately. Making 
use of Eq. lj2U|) we get: 

5W {ah) (t) = - — / dxdydz GST*, (27) 

7T 2 Jo 

where 



x y 

C = cosuj(t — x — y — z) cosuj ( t — — — — 



(x z \ / y z 
t ] cos u> 1 1 
2 2/ V 2 2 



§ = 3sin 2 ^ _ sin 2 *v _ sin 2 tl - 4 rsin ^ 
2 2 2 2 



x sin -jj- sm — (x sm xr x + y sin 



: sin 



D is the product of three diffusons in Eq. (|21ll evaluated 
at r\ = 0, and we introduced $ x = y — z, {ty — —x — z, 
and d z — x + y. 

The long-time behavior of Eq. I|27|l is determined 
by the vicinities of the no-dephasing points |1(J) where 
each of the three diffusons entering D is equal to 1. 
An analogous situation arises in the calculation of the 



one-loop quantum correction for the periodically driven 
orthogonal matrices @|, which is dominated by the no- 
dephasing points of a single dynamic cooperon. In 
the present case, the no-dephasing points are given by 
(x, y, z) = (x, 2-nm/uj — x, 2im/uj — x) with arbitrary x 
and integer m and n. 

In the limit t ^> (w _1 , the no-dephasing points 
with different m and n do not overlap and the triple 
integral in Eq. lt2"?T) can be evaluated as 



dx dy dz 



dxj:J 



dSy dSz, 



(28) 



where we introduced y = 2 r Kmjio — x + 5y and z = 
2im/uj — x + Sz. At the no-dephasing point the factor C 
is nonzero whereas the factor 8 vanishes and should be 
expanded in the deviations Sy and Sz: 

6 = cos 2 ujt cos 2 u(t + x), (29) 

§ = ^(3(Sy~6z) 2 -6z 2 -Sy 2 ) 
Tuj 4 

H —SySz(Sy — Sz) [(x + z)Sy — (x + y)Sz] . (30) 

Though the last term of Eq. (|3U|) is proportional to 
the fourth power of Sy and Sz, their smallness is com- 
pensated by an extra factor x,y, z ~ t. In the limit 
t ^> (uj^ 1 , r _1 ) we can integrate near the no-dephasing 
points in the Gaussian approximation retaining only 
quadratic in the deviations terms in In D: 

T> = cxp j-^- [x(8y - Sz) 2 + ySz 2 + zSy 2 ] J . (31) 

The weight (|31(l determines the correlators: 

M _ ( (SySy) (SySz)\ 
\(SySz) (SzSz)) 



1 



x + y 



Flo 2 xy + yz + zx 



(32) 



Substituting Eqs. into Eq. gZJ and inte- 

grating over Sy and Sz one gets 



SW {ab \t) 



■ COS 2 Lot 



x dx 



mn 



27iVdetAf (§), (33) 



where we replaced cos 2 u>(t + x) by its average value 1/2. 
The average (§) is calculated with the help of the Wick's 
theorem using the pair correlators l|32l) : 



(§) = 



"ixyz 



2T(xy + yz + zx) 2 



(34) 
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Finally, since the summand in Eq. (|33[) is a smooth func- 
tion of to and n it is possible to pass from summation 
over 777 and 77 back to integration over y and z: 



E 

inn 



(- 

V27T 



dy dz. 



(35) 



As a result we obtain 



5W {ab \t) = - 



ZAlo 2 



4?T 3 



■ cos cot 

x+y+z<t 



xyz dx dy dz 



(xy + yz + zx) 5 / 2 
This integral is equal to (2w/27)t and we get 

Au 2 t 



8W {ab) {t) = 



18tt 2 



■ COS LOt. 



(36) 



(37) 



The contribution of the diagram (c), 8W(°\ can be 
calculated analogously. Due to the same structure of 
the diffusons, its no-dephasing points coincide with the 
no-dephasing points for 5W^ ab \ Instead of Eq. Q33JI one 
has now: 



8W {c) {t) = 



TAlo 2 d 
Air 2 dt 



dx 



mn 



27rVdetM(S'), (38) 



where 



(§') = 1 - Flo 2 {(5y - Sz) [(x + z)5y - (x + y)5z\) 



yz 



xy + yz + zx 



(39) 



Passing from summation to integration according to 
Eq. (|35|l and utilizing the symmetry properties of the 
integrand we obtain: 

{B *) (t) = _£4| f" t * + « (40) 
; " ' sjxy + yz + zx y ' 



24tt 3 dt J 



The integral is equal to (ir/6)t 2 yielding 



SW (c) (t) 



Ato 2 t 



(41) 



Note a peculiar property of Eqs. I|37|l and l|41|) : 

SW^^t) oc t(dif/dt) 2 and vanishes at the turning 
points of the perturbation, whereas 5W^(t) is always 
positive, even when dip/dt = 0. This means that they 
describe different mechanisms of absorption, with dif- 
ferent memories on the past. 

Combining Eqs. (|3"7|l and i|4T|) we get the total two- 
loop correction to the quasiclassical absorption rate in 
the harmonic case: 



SW(t) 



Auj 2 t 
' 72tt 2 



[4 cos 2 tot + l] 



(42) 



valid at t > (a; -1 , T" 1 ). 

The time-averaged correction grows linearly with the 
duration of the perturbation: 



SW(t) = - 



AloH 
24tt 2 ' 



(43) 



Remarkably, Eq. I|43|l holds not only for a harmonic 
perturbation but for a?7 arbitrary periodic perturbation 
with the period 2tt/lo. Formally this follows from the 
fact that the level sensitivity T to the external perturba- 
tion drops from Eq. (|43|l . Then, according to Eq. (|35|l . 
the factor to 2 in Eq. I|43|) measures the inverse time sep- 
aration between the no-dephasing points which is the 
same for all periodic perturbations of a given period. 

6. Dynamic vs. Anderson localization. It is 
useful to compare the two-loop result 1)43(1 for a har- 
monic perturbation with the analogous one-loop expres- 
sion for the GOE obtained in Ref. 0]: 



SW(t) 

Wn 



-, GOE, 



7rf 

24L 



(44) 



, GUE, 



where Wq = ttTlo 2 /2A is the period-averaged absorp- 
tion rate, and £* = ?T 3 r/2A 2 is the localization time. 

In Ref. 0] we pointed out that the weak dynamic 
localization correction to the energy absorption rate of 
a periodically driven GOE has the same square-root be- 
havior as the weak Anderson localization correction to 
the conductivity of a quasi-one-dimensional (ID) disor- 
dered wire. Now we see that the same is true for the 
case of the GUE as well: in both cases the correction is 
linear in time and dephasing time, respectively. There- 
fore it is tempting to suggest that this analogy is not a 
coincidence but has its roots in equivalence between the 
dynamic localization for the RMT driven by a harmonic 
perturbation and ID Anderson localization. 

Such an equivalence is known for the case of kicked 
quantum rotor (KQR): in the long time limit, the KQR 
problem can be mapped m onto the ID er-model. On 
the other hand, the problems of the <5-kicked KQR and 
of the periodically driven RMT are, to some extent, 
complementary. Both of them can be mapped on a 
tight-binding ID model, but with very different struc- 
ture of couplings between the sites and auxiliary or- 
bitals 0]. In particular, the "kicked RMT" model with 
ip(t) being a periodic S- function does not exhibit dy- 
namic localization whatsoever 0]. 

In order to check the assumption about the equiv- 
alence of the driven RMT to the quasi- ID disordered 
wire we use the simple relationship between the time- 
dependent energy absorption rate W(t) in the dynamic 
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problem and the frequency-dependent diffusion coeffi- 
cient D(ui) in the Anderson model |12j : 



W(t) 

w 



-t-oo 



2tt 



-iuj ■ 



' DM 
D ' 



(45) 



where Wo and D are the classical period-averaged ab- 
sorption rate and diffusion coefficient. D(ui) is known 
from the theory of weak Anderson localization: 



1 



SD(lu) 
Do 



y/-iuti 
f 



6iuit\ oc 



GOE, 
GUE. 



(46) 



Here t\ oc = (27Wi) 2 Z?o, and v\ is the ID density of 
states. ThenEqs. (gBJ), (gSJl give two expressions similar 
to Eq. (|44|l with only one fitting parameter t*/t\ oc . One 
can easily see that with the choice t*/t\ oc = ir/A both 
numerical coefficients match exactly. 

We believe that there are deep reasons for this 
coincidence and make a conjecture that the (period- 
averaged) dynamics of the harmonically-driven RMT at 
time scales t ^> (a; -1 , T _1 ) is equivalent to the den- 
sity propagation in a quasi-fD disordered wire. Em- 
ploying this equivalence, we can easily calculate the en- 
ergy absorption rate in the regime of well developed dy- 
namic localization at t 3> using the Mott-Berezinsky 
asymptotics of the AC conductivity, o~(uj) cx ui 2 ln 2 (f /lo) 
[1311111. Substituting D(uj) cx a(w) into Eq. gHJ) we find 
that in the localized regime W(t) decays as 



W{t) cx 



lnt 



t > u 



(47) 



This dependence is not directly related to the spa- 
tial dependence of the localized wave functions which 
is exponential in the Anderson model. It can be 
seen if one considers the density-density correlator 
[disorder-averaged product of the retarded and ad- 
vanced Green's functions G R (x, x' , e + uj)G a (x', x, e)] 
whose Fourier transform can be conveniently repre- 
sented as 2-kvx A(k, uj)/(— iuj). According to Gorkov's 
criterion of localization [TJ], A(k,0) is finite and its 
Fourier transform determines the spatial decay of lo- 
calized wavefunctions. On the other hand, D(oj) can be 
extracted from the density-density correlator as 



iuj d 2 
Tdk? 



A(k,u) 



k=0 



(48) 



and, according to our conjecture, should be substituted 
in Eq. I|45|l to give the absorption rate. Thus, instead of 
A(k,oj — 0), usually studied in the Anderson localiza- 
tion problem, W(t) is determined by the to dependence 



of d 2 A{k,uf) / dk 2 at k = 0, which to the best of our 
knowledge evaded investigation in the framework of the 
quasi- ID nonlinear sigma model. 

7. Conclusion. We derived the general expression 
for the lowest order (two-loop) interference correction 
to the energy absorption rate of a parametrically-driven 
GUE. If an external perturbation grows linearly with 
time the first correction vanishes. For a periodic pertur- 
bation the averaged correction SW(t) oc t. We make a 
conjecture that the dynamics of the harmonically-driven 
RMT at the time scales t 3> 1/T is equivalent to 
the ID Anderson model. Based on this equivalence we 
predict that in the regime of strong dynamic localization 
W(t) cx \n{t)/t 2 . 
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